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Abstract —In this paper, a new mnlti-cell MMSE precoder is 
proposed for massive MIMO systems. We consider a multi-cell 
network where each cell has K users and B orthogonal pilot 
sequences are available, with B = PK and /3 > 1 being the pilot 
reuse factor over the network. In comparison with conventional 
single-cell precoding which only uses the K intra-cell channel 
estimates, the proposed multi-cell MMSE precoder utilizes all B 
channel directions that can be estimated locally at a base station, 
so that the transmission is designed spatially to suppress both 
parts of the inter-cell and intra-cell interference. To evaluate 
the performance, a large-scale approximation of the downlink 
SINK for the proposed multi-cell MMSE precoder is derived 
and the approximation is tight in the large-system limit. Power 
control for the pilot and payload, imperfect channel estimation 
and arbitrary pilot allocation are accounted for in our precoder. 
Numerical results show that the proposed multi-cell MMSE 
precoder achieves a significant sum spectral efficiency gain over 
the classical single-cell MMSE precoder and the gain increases as 
K or p grows. Compared with the recent M-ZF precoder, whose 
performance degrades drastically for a large K, our M-MMSE 
can always guarantee a high and stable performance. Moreover, 
the large-scale approximation is easy to compute and shown to 
be accurate even for small system dimensions. 

1. Introduction 

Massive multiple-input-multiple-ouput (MIMO) is an excit¬ 
ing new multi-user MIMO technology that has gained traction 
in recently years Cl, d. The idea is to employ an array 
comprising hundreds of antennas at each base station (BS) 
and serve tens of users simultaneously. The system spectral 
efficiency (SE) can be improved by an order of magnitude 
without consuming extra spectral resource ci-ii. By using 
simple coherent linear precoders such as zero forcing, the 
intra-cell interference and receiver noise can be averaged out 
in the limit of an infinite number of BS antennas, and the only 
remaining performance bottlenecks are pilot contamination 
and distortion noise from hardware impairments d, 0 . The 
uplink and downlink transmit powers can also be reduced by 
an order of magnitude since a comparable array gain can be 
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harnessed by the phase-coherent processing 16]. Eurthermore, 
in time division duplex (TDD) mode, the channel estimation 
overhead scales linearly with the number of users, instead of 
the number of BS antennas, which allows for adding extra 
antennas without affecting the estimation overhead C). These 
features make massive MIMO one of the key technologies for 
the next generation wireless communication systems. 

In the downlink, the most common linear precoding 
schemes are matched filtering (ME), zero forcing (ZE) and 
minimum mean square error (MMSE)Q We consider a multi¬ 
cell network with B orthogonal uplink pilot sequences and K 
users in each cell, where B = PK and /3 > 1 is called the 
pilot reuse factor (i.e., only 1//3 of the cells use the same set 
of pilots). In conventional massive MIMO systems, each BS 
estimates the K intra-cell channels by listening to the pilot 
signalling from its own cell, and then constructs user-specific 
precoders based on the channel estimates to limit intra-cell 
interference ID, 0, Hi- When /3 > 1, however, the BS can 
locally estimate B, instead of K, channel directions by also 
listening to the pilot signalling from all cells. Thereby, the BS 
is able to select its precoders to actively suppress also the inter¬ 
cell interference, since its K serving users only occupy K out 
of the B channel directions. A downlink multi-cell precoding 
scheme is proposed in ifTO] , which brings a notable gain over 
single-cell processing. However, this scheme does not account 
for optimized pilot allocation which, as shown in ID, is an 
important way to suppress pilot contamination and achieve 
high system SE in massive MIMO deployments. Moreover, 
no closed-form performance analysis is provided in ifTo] . 

In this paper, we propose a new multi-cell MMSE pre¬ 
coder and derive a large-scale approximation of the downlink 
SINK, which is tight the large-system limit. Power control 
for the pilot and payload, imperfect channel estimation and 
arbitrary pilot allocation are all accounted for in our scheme. 
By utilizing all the channel directions that can be locally 
estimated at a BS, the proposed multi-cell MMSE precoder 
can not only suppress intra-cell interference, but also actively 
reduce the interference caused to other cells. Numerical results 
are provided to show that a significant sum SE gain can 

*A special case of the MMSE precoder is the regularized ZF precoder, 
which is obtained when all the users in a cell have equal pathlosses. 


be obtained by the proposed precoder over the single-cell 
MMSE precoder and the gain becomes more signihcant as 
K or /3 increase. Compared to the multi-cell ZF precoder 
from a, whose performance degrades drastically for larger 
K or P, our M-MMSE precoder provides a high and stable 
performance even for very large K. Furthermore, the large- 
scale approximation of the SINR is easy to compute and shown 
to be a very accurate approximation even for small system 
dimensions. 

Notations: The trace, conjugate, transpose, Hermitian trans¬ 
pose and matrix inverse operators are denoted by tr(-), (•)*, 
(■)^’ respectively. 

II. System Model and Precoder Design 

We consider a massive MIMO cellular network operating 
according to a synchronous TDD protocol. Each cell is as¬ 
signed an index in the cell set £, and the cardinality |£| is 
the number of cells. The BS in each cell is equipped with M 
antenna elements and serves K single-antenna users within 
each coherence block. Assume that this time-frequency block 
consists of Tc seconds and Wc Hz, such that Tc is smaller 
than the coherence time of all users and Wc is smaller than 
the coherence bandwidth of all users. This leaves room for 
S = TcX Wc transmission symbols per block, and the channels 
of all users remain constant within each block. Let hjik denote 
the channel response from user k in cell I to BS j within a 
block and assume that it is a realization from a zero-mean 
circularly symmetric complex Gaussian distribution: 

hjifc ^ CA/’(0, dj (zjfc) Im) . (1) 

The vector zik G is the geographical position of user k in 
cell I and dj (z) accounts for the channel attenuation (e.g., path 
loss and shadowing) related to some user position z. Since the 
user positions change relatively slowly, dj (zik) is assumed to 
be known at BS j for all I and all k. 

In the TDD protocol, the channels are estimated from 
uplink pilot signaling and utilized for downlink precoding by 
exploiting channel reciprocity. Thus, the channel estimation is 
first discussed to lay a foundation for the precoder design. 


assume that all pilot sequences originate from a predefined 
orthogonal pilot book, defined as V = {vi,..., vs}, where 




B, bi = &2, 

0, bi ^ b2, 


(3) 


and let ilk G {1,... ,B} denote the index of the pilot sequence 
used by user k in cell 1. Arbitrary pilot allocation is supported 
and the relation between B and K is denoted hy B — j3K, 
where /3 > 1 is called the pilot reuse factor. If the pilots are 
allocated wisely between the cells, a larger /3 brings higher 
estimation quality and a lower level of pilot contamination. 

Based on the received pilot signal in dal), the MMSE 
estimate of the uplink channel h.jik is 

= V^kd, {zik) Y, (vl>*) V*^, (4) 


where is the covariance matrix of the vectorized received 
signal vec (Yj) and is given by 

K 

Pimdj (Zim) -b Cr^Is. (5) 

r6£m=l 


Then according to the orthogonality principle of MMSE 
estimation, the covariance matrix of the estimation error 
hjik = h^ik - hjtk is given by 

Cjik = IE I 

= dj (zik) (1 - pikdj (zik) Im 

= dj{zik){l-pikdj{zik)aji^^B)lM, ( 6 ) 

where and can be simplified as 


J'^lk V ^ ^ \ H I 2 * ' ^ 

As pointed out in H, the term Yj ^ v*^ in dUl 

depends only on which pilot sequence that user k in cell I 
uses. Consequently, users who use the same pilot sequence 
have parallel estimated channels at each BS, while only the 
amplitudes are different. To show this explicitly, define the 
M X B matrix 


A. Uplink Channel Estimation 

In the channel estimation phase, the collective received 
signal at BS j is denoted as Yj G where B is the 

length of a pilot sequence (and thus the number of orthogonal 
pilot sequences available in the network). Then can be 
expressed as 


K 

= E E , ( 2 ) 

fc=l 

where h.jik is the channel defined in ([T]i and pik is the power 
control coefficient for the pilot of user k in cell 1. The matrix 
Nj contains independent elements which follow a complex 
Gaussian distribution with zero mean and variance We 


Hv., = = Y, {^*) [vt,...,v^], (8) 

so that the MMSE estimate in (|4|i can be expressed as 

^jlk — s/Pikdj (z^/i;) Hy , (9) 


where denotes the ith column of the identity matrix I^. 
The fact that users with the same pilot have parallel estimated 
channels is used herein to derive a new SE expression. 

Notice that the estimated channel hjik is also a zero-mean 
Gaussian vector, and its covariance matrix ^jik is 

^jik = dj (z/fe) Im — Cjik = pikd^j (z/^) aji^^Blj^. (10) 

Define the covariance matrix of hyji as Then accord¬ 

ing to ® and (fTOl) . = ctjiBlM- 





Vjk — 


(15) 
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E{h} < 
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B. Downlink Multi-cell MMSE Precoder 


downlink SE 


During the downlink payload transmission, the collective 
received signal at user k in cell j can be expressed as 


Rjk — 



log 2 (1 + T]jk ) 


(14) 


K 

Vjk Veim^lrnSlm + Tljk, (11) 

TTi — l. 


where wim S is the precoder associated with user m 

in cell I, sim CA/^(0,1) is the corresponding payload data 
symbol, qi^ is the downlink transmit power coefficient, and 
rijk ^ CAf (0,1) is additive white Gaussian noise (AWGN). 

In ifTTl . we propose a multi-cell MMSE (M-MMSE) detec¬ 
tor for the uplink of massive MIMO. In contrast to conven¬ 
tional detectors that only utilizes the K intra-cell channels, our 
proposed M-MSME detector utilizes all B estimated channel 
directions in the matrix Hy j from (H). Since the K users only 
occupy K out of the B channel directions, our M-MMSE 
detector can suppress both intra and inter-cell interference, 
whereas the conventional detectors can only suppress intra¬ 
cell interference. The proposed M-MMSE detector maximizes 
the uplink SINR under general conditions, and notable SE 
gains are achieved over conventional single-cell MMSE and 
the multi-cell ZE. The M-MMSE detector from UTl is given 
by 


gM-MMSE ^ + (^2 ^ 

( 12 ) 

K 

where = X] BkPikd'j and Tik is the 

l^L k—1 

Uplink transmit payload power from user k in cell 1 . Moreover, 

K 

Pj = Bkdj (zik) (l-pikdj (Zifc) aji.^B), where 

lec k=i 

is dehned in (|7]i. 

Recently, a established an uplink-downlink duality for 
massive MIMO systems which proves that the uplink SEs can 
also be achieved in the downlink if each downlink precoder 
is a scaled version of the corresponding uplink detector. As 
the M-MMSE detector is the state-of-the-art uplink method, 
we propose to apply the same methodology for downlink 
precoding. The downlink M-MMSE precoder is constructed 
as 


M-MMSE 


M-MMSE 

_ 

\/^jk 


(13) 


where normalizes the average trans¬ 

mit power for the user k in cell j to ^{W-s/Qh^WjkSimW^} = 
gim- Since there are no downlink pilots in our TDD protocol, 
the users do not know the current channels but can only 
learn the average channel gain, y^p^Eh{hj^^Wjfc}, and the 
total interference variance. Therefore, the ergodic achievable 


can be obtained at user k in cell I as in Q, lU, where (1 — ^) 
is the pilot overhead and rjjk is given at top of this page. 

The downlink SINR in ( fTSl l holds for any linear precoding 
scheme. This SE is achieved by treating E|h}{hj^^Wjfc} as 
the true channel in the receiver, and treating interference and 
channel variations as worst-case Gaussian noise. Thus, Rjk is 
a lower bound on the downlink ergodic capacity. 

By utilizing all the available estimated directions, the M- 
MMSE precoder suppresses both intra-cell interference and 
the interference caused to other cells, and thus a higher SINR 
can be expected by our precoder than conventional single¬ 
cell precoders, at least for an appropriate power control (as 
prescribed by the uplink-downlink duality concept in H). In 
the next section, a large-scale approximation of the downlink 
SINR in (flST l is derived. Note that in ifTOll . the authors also 
propose a multi-cell MMSE precoder, which brings a notable 
gain over single-cell processing. But it does not accounted 
for arbitrary or optimized pilot allocation which, as shown 
in a, is an important way to suppress pilot contamination 
and achieve high system SE in massive MIMO. Moreover, no 
closed-form performance analysis is provided in nni. 

III. Asymptotic Analysis 

In this section, performance analysis is conducted for the 
proposed M-MMSE precoder. The large-system limit is con¬ 
sidered, where M and K go to inhnity while keeping the 
ratio K/M hnite. In what follows, the notation M ^ oo 
refers to AT, M —?> oo such that lim sup^AT/M < oo and 
liminfMAT/M > oil Since B scales with K for a hxed /3, 
limsupjv^S/M < oo and liminfMf?/Af > 0 also hold for B. 
The results should be understood in the way that, for each set 
of system dimension parameters M, K and B, we provide an 
approximative expression for the SINR, and the expression is 
tight as M, K and B grow large. In what follows, the notation 
- > denotes the convergence of a deterministic sequence. 

M^oo 

Before we continue with our performance analysis, two use¬ 
ful results from large random matrix theory are first recalled. 
All vectors and matrices should be understood as sequences 
of vectors and matrices of growing dimensions. 

A. Useful theorems 

Theorem 1 (Theorem 1 in ifT^ I: Let D s ^mxm 
terministic and H G £_mxb |,g with independent 

^The limit superior of a sequence Xn is defined by limsup^a:^ = 
lim (sup {xm : rn ^ n}); the limit inferior is defined as liminf^x^ = 

n—>-oo 

lim (inf {xm '■ ttl ^ n}). 

n—^oo 
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column vectors ^ CAf (O, Assume that D and the Theorem 3 For the downlink MMSE precoder in ( fTSl) . we 

matrices R;, (b = 1, have uniformly bounded spectral have rjjk — fjjk ->■ 0, where fjjk is given at top of this 


norms (with respect to M). Then, for any p > 0, 


1 


rtr 


D (hh" + pIm) 


-1 
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where T (p) G C 
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(16) 


1 / ' \ 


^Ijkm = J 

(17) 

Pjkdl (Zjfc) ^li. 


2 + 'yiiji^'dijk 
(l + 'Ylijk'dljk) 


(26) 


(27) 


and the elements of 6 (p) = [5i (p), Sb (p)]^ are dehned 
as Sb (p) = limt^.oo (p) ,b = I,B, where 

/ 


(P) = ^tr 
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Rb — y 

\ M ^ 


R, 


V 




1 + (p) 


pi 


■N 


for f = 1, 2,..., with initial values = 1/p for all b. 


(18) 


Theorem 2 (see ifT^ l Let 0 G i^^mxm Hermitian non¬ 
negative dehnite with uniformly bounded spectral norm (with 
respect to M). Under the same conditions for D and H as in 
Theorem [T] 

y tr (b (hH^ + pIm) © (hH^ + pIm) 

-ytr(DT'(p))-^0 (19) 

Ad M—¥00 

where T' (p) G is dehned as 

T' (p) = T (p) 0T (p) + T (p) y y (p). 

( 20 ) 

T (p) and S (p) are given by Theorem [T] and S' (p) = 
(p)]^ is calculated as 

5'(p) = (Ib-J(p))-'v(p) (21) 


where J (p) and v (p) are dehned as 


[J(p)]m = 


^tr(R,T(p)R;T(p)) 

M{l+Si{p)f 


,l<b,l< B 


[v (p)], = y tr (RbT (p) 0T (p)), 1 < 5 < 5. 


( 22 ) 

(23) 


B. A Large-Scale Approximation of the SINK 

In what follows, we derive the large-scale approximation 
fjjk of Pjk with the M-MMSE precoder such that 


Vjk - Vjk -0. (24) 

M—¥00 


"dijk = ^tr , (28) 

y-fem = , (29) 

■d'lm = > (30) 

where 

1) T; = T; (a) and S^a) = [(5i,are given by 

Theorem [U for a = and R;, = 7/b€*v,z&- 

2) ^'im = T^'im (a) and S' (a) = are given by 

Theorem |2] for a = 0 = and Rf, = 

Vlb^Vjb- 

3) T;'„ = t"^ (a) and^^' (a) = ...,6'gf are given by 

Theorem|2]for a = © = Im, and R^ = 7ib$v,zb- 

Proof: See the Appendix B and C of ifTJl . ■ 

By utilizing Theorem [3l the ergodic SE Rjk in (fT4l) . after 
dropping the prelog factor, converges to Rjk = log 2 (l+p?fc) in 
the large-system limit. Therefore, a large-scale approximation 
of the downlink ergodic SE is provided by (1 — ^)Rjk- This 
approximation is easy to compute and, as shown in SectionITVl 
is very accurate even at small system dimensions. 

IV. Simulation Results 

In this section, we illustrate the accuracy and usefulness 
of the analytical contributions for a symmetric hexagonal 
network topology. We consider the classic 19-cell-wrap-around 
structure to avoid edge effects and guarantee the same SE in 
all cells a, 113 . Each hexagonal cell has a radius of r = 500 
meters. To achieve a symmetric pilot allocation network, the 
non-universal reuse factor can be /3 G {1, 3,4, 7, • • • }, and the 
pilots are then allocated randomly within each cell. 

User locations are generated independently and randomly in 
the cells by following uniform distributions, but the distance to 
the serving BS is at least 0.14r. Eor each user location z G K^, 
a classic pathloss model is considered, where the variance of 
channel attenuation is dj (z) = H^re G is the 

location of BS in cell j, k is the pathloss exponent, and H-H 
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Fig. 1. Achievable sum SE as a function of M, for /3 € {1,3,4, 7} and 
K = 10. 

denotes the Euclidean norm. C > 0 is independent shadow 
fading with lOlog^Q (C) ^ A/’(0,cr^j). We assume k = 3.7, 
= 5, and coherence block length S = 5000 

Statistical channel inversion power control is applied to 
both uplink pilot transmission and the “virtual” uplink pay- 
load power used in the M-MMSE precoder construction, 
i.e., pik = Tik = a7(i7J) ISI- Thus the average uplink 
SNR per antenna and user at its serving BS is constant: 
^{pik ||hnfc||^ }/(Mcr^) = pja^. This is a simple but effective 
policy to avoid near-far effects. Eor downlink payload data 
transmission, equal power allocation is used for each user, i.e., 
gik = Pmax- In the simulations, p/cr^ is set to OdB to allow 
for decent uplink channel estimation accuracy, and Pmax is 
selected to make the cell edge SNR (without shadowing) equal 
to — 3dB in the downlink. 

To verify the accuracy of the large-scale approximations, 
10000 independent Monte-Carlo channel realizations for small 
scale fading are generated to numerically calculate the down¬ 
link achievable SE in (fl4ll . The numerical result and its 
approximation from Theorem |3] are shown in Eig. □ As 
seen from the figure, the achievable sum SE increases mono- 
tonically with (3 for the considered range of values. This 
is due to the following two properties. Eirstly, a larger /3 
results in a lower level of pilot contamination, contributes to a 
higher channel estimation accuracy, and thereby increases the 
achievable SE. Secondly, a larger /3 makes it possible for the 
BS to estimate more channel directions to users in other cells 
and design the precoders to reduce the interference caused 
to these other-cell users. However, it should be notice that a 
larger (3 results in a larger pre-log loss, and thus the achievable 
SE may first increase and then decrease as (3 grows when S 
is small. Eig. [T] also shows that the numerical results and the 
large-scale approximations match very well, even for relatively 
small M and small K. 

To show the advantages of our M-MMSE precoder over 
other schemes, simulation results for the matched filter (ME) 

2 This coherence block can, for example, have the dimensions of Tc = 5 ms 
and Wc = 100 kHz. 


Fig. 2. Achievable sum SE of M-MMSE, M-ZF, S-MMSE and M-MF for 
/3 = 4. 

from m , the multi-cell ZE (M-ZE) precoding from a, and the 
S-MMSE precoding from jS] are provided for comparison in 
Eig. 12] The same power allocation and normalization from (fOb 
are applied for all precoders. Note that M — I3K > 0 is 
needed for the M-ZE scheme, thus the minimum value of M 
for the M-ZE is f3K + 1. We provide the results for /3 = 4 
in Eig. 121 As seen from the figure, the ME scheme achieves 
the lowest performance since it does not actively suppress 
any interference. Compared with the S-MMSE, our M-MMSE 
achieves a notable SE gain and the gain increases as K grows. 
Eor M = 200, the SE of M-MMSE is 10% higher than that of 
S-MMSE with K = 10, and the gain increases to 20% with 
K = 30. The advantage of the M-MMSE over the M-ZE is 
relatively small when K = 10, but it becomes significant as 
K grows. The performance gains over the S-MMSE and the 
M-ZE increase as (3 grows as well, and the simulations results 
are provided in ifTSl due to the space limitation. Moreover, the 
M-ZE can sometimes achieve very low SE for small M, while 
our M-MMSE always seems to achieve good performance. In 
addition, the computational complexity of the M-MMSE is 
the same as for the M-ZE. Thus in general, our M-MMSE 
precoder is the better choice if high system SE is desirable. 

To show the ability of our M-MMSE precoder to handle 
different number of users, the SE is shown in Eig. [3] as a 
function of K for M = 200. Eor each K, /3 G {1, 3,4, 7} are 
considered, and for each set of K and f3, we use the large- 
scale approximation given in Theorem [3] to obtain the sum SE 
for the M-MMSE precoder. Eig. [3] shows the SE with the [3 £ 
{1,3,4, 7} that that gives the highest SE for each K. The result 
for M-ZE is provided for comparison. As seen from the figure, 
the optimal value of (3 changes for different K. Taking the M- 
MMSE precoder as an example, /3 = 7 gives the highest SE 
for K < 20, then /3 = 3 is desirable for a wide range of values 
of K, and finally the best choice is /3 = 1 for K > 100 in 
order to reduce the pilot overhead. The maximal SE is achieved 
by /3 = 3 for both precoders, but the corresponding K and 
SE are higher for M-MMSE. It is notable that our M-MMSE 
can achieve high SE even for hundreds of users, while the 
performance of the M-ZE degrades drastically when K > 70, 
























Fig. 3. Achievable sum SE as a function of K with M = 200. 


mainly since the array gain M — j3K is much smaller than M 
in this range. Therefore, our M-MMSE can guarantee a much 
higher performance than the M-ZF in a crowded network. 

Finally we consider the rehned pilot reuse scheme proposed 
in Cl where pilot reuse is only applied at the cell edges. This 
reuse policy selects B = Kf + {1 — /3/)/3) pilot sequences, 
and assign the same PfK pilot sequences to the PfK users 
that are closest to their serving BSs. The other K{1 — I3f )fi 
pilots are assigned to the remaining users according to the 
same symmetric reuse policy as before. The rehned scheme 
reduces to the symmetric reuse policy for /?/ = 0, while /3/ = 
1 means that the same pilots are used in all cells. The idea of 
this scheme is to mitigate pilot contamination while limiting 
the pilot overhead. We apply this scheme to our M-MMSE 
precoder, and show the performance in Fig. |4] for M = 200 
and different K. For each K and /?/, we choose the highest 
SE from what we obtained from /3 G {1,3,4, 7}. As seen 
from the hgure, /3/ = 0 is preferred when K is small, which 
means that pilot reuse should also be applied for everyone. 
However, ,8/ > 0 is benehcial for larger K values since there 
will likely be cell center users that are not sensitive to pilot 
contamination, which can be exploited as in Cl to reduce 
the pilot overhead. However, if we look at the highest SE for 
each Pf among the considered K and /3, the increase of /?/ is 
not as obvious as it is for M-ZF shown in M . thanks to the 
strong ability of our M-MMSE to actively suppress inter-cell 
interference. 

V. Conclusions 

In this paper, a multi-cell MMSE precoder is proposed 
based on uplink-downlink duality, and an approximative SINR 
expression is derived which is tight in the large-system limit. 
Compared with conventional single-cell precoders, that only 
make use of the estimated channel directions within the 
serving cell, our multi-cell MMSE scheme utilizes all channel 
directions that can be estimated at the BS to also mitigate inter¬ 
cell interference. Numerical results show that the proposed 
precoder brings a very promising sum SE gain over the single¬ 
cell MMSE and the gain becomes signihcant as AT or /3 grows. 
Compared to the M-ZF, our M-MMSE precoder can maintain 



Fig. 4. Achievable sum SE as a function of (if with M = 200. 

high and stable SE for very large K. Furthermore, the SINR 
approximation is easy to compute and shown to be accurate 
even for small system dimensions. 
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